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Aims. The goal of this study is to perform numerical tests of Zeeman Doppler Imaging (ZDI) to asses whether correct reconstruction 
of magnetic fields is at all possible without taking temperature into account for stars in which magnetic and temperature inhomo- 
geneities are spatially correlated. 

Methods. We used a modern ZDI code employing a physically realistic treatment of the polarized radiative transfer in all four Stokes 
parameters. We generated artificial observations of isolated magnetic spots and of magnetic features coinciding with cool temperature 

O spots and then reconstructed magnetic and temperature distributions from these data. 
Results. Using Stokes / and V for simultaneous magnetic and temperature mapping for the star with a homogeneous temperature dis- 
tribution yields magnetic field strengths underestimated by typically 10-15% relative to their true values. When temperature is kept 
constant and Stokes / is not used for magnetic mapping, the underestimation is 30-60%. At the same time, the strength of magnetic 
field inside cool spots is underestimated by as much as 80-95% and the spot geometry is also poorly reconstructed when temperature 
■ variations are ignored. On the other hand, the inversion quality is greatly improved when temperature variations are accounted for 

in magnetic mapping. The field strength is underestimated by 40-70% for the radial and azimuthal spots and by 70-80% for the 
meridional spots. Inversions still suffer from significant crosstalk between radial and meridional fields at low latitudes. The azimuthal 
field component proves to be most robust since it suffers the least from crosstalk. When using all four Stokes parameters crosstalk is 
removed. In that case, the reconstructed field strength inside cool spots is underestimated by 30-40% but the spot geometry can be 
recovered very accurately compared to the experiments with circular polarization alone. 

Conclusions. Reliable magnetic field reconstruction for a star with high-contrast temperature spots is essentially impossible if tem- 
perature inhomogeneities are ignored. A physically realistic line profile modeling method, which simultaneously accounts for both 
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types of inhomogeneities, is required for meaningful ZDI of cool active stars. 
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^\ , 1 . Introduction a visible pole will only produce a distortion close to the middle 

OO ■ of the line profile. 

q ■ Stellar magnetic fields vary widely in strength and configuration. If high-resolution circularly polarized (Stokes V) observa- 
~- 1 There are many types of magnetic field generation mechanisms, tlons are included in DI, a tool for magnetic field reconstruc- 
O ; but not even the magnetic field generation in the Sun is com- tlon < called Zeeman Doppler Imaging (ZDI) < |Semel| | 1989|) , is 
rl ■ pletely understood. To be able to understand the different types obtained. It uses the same basic principles as conventional DI, 
of magnetic fields and their origin we first need to know how the but taes t0 solve a far more challenging problem of interpret- 
fields are structured on the stellar surfaces and how their config- m 8 polarization signatures inside spectral lines in terms of the 
urations evolve with time surface distribution of vector magnetic fields. 

N owadays ZDI is a widely used t echnique dDonati et al.l 

Doppler Imaging (DI) has long been used to reconstruct 2003; Petit et al. 2004; Catala et al. 2007), applied to many dif- 
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' 2D maps of the distribut ion of temperature or chemical abun- ferent types of hot and cool magnetic stars. This widespread 

dance on stellar surfaces (IVogt et alll 9871: iPiskunov et al.|[T990l usage requires testing validity of the method. Several studies 

Koch ukhov et al. 2004). This technique uses rotationally modu- demonstrated t hat ZDI is u s eful, a lthough it has its limitations, 

lated unpolarized (Stokes /) line profiles at different rotational For example, Bro wn et"a"fl dl99ll) found that ZDI using only 

phases to reconstruct a map of inhomogeneities on the stellar Stokes V could reproduce a spotted magnetic field structure 

surface. Temperature spots on cool active stars will show up in more accurately than a large dipole-like structure. They did how- 

the line profiles as intensity variations, or bumps. The position ever chose to dissociate temperature spots and magnetic spots 

of the temperature spot on the star will determine the position of by reconstructing a brightness map for a star with temperature 

the bump within the line profile. Reversely, the position of the spots only, and, conversely, by only reconstructing the magnetic 

bump can be used to i nfer the longitude of the temperature spot field for a star with magnetic spots but without any tempera- 

dVogt & Penrodll9"83h . By analyzing line profiles at different ro- ture inhomogeneities. They also used an analytical line forma- 

tational phases, the latitude of the temperature spot can also be tion model assuming a plane parallel Milne-Eddington model 

determined. If the spot lies close to the rotational equator, the atmosphere and applied inversions to simulated observations of 

bump will show up in the far blue wing of the profile and then a single magnetically sensitive spectral line. Each magnetic map 

move across the line. On the other hand, a temperature spot near was reconstructed under the assumption of a radial field. 
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Further tests carried out bv lDonati & Brownl (1 19971) showed 
that ZDI can, to some extent, reproduce the field orientation in 
a spot. The method proved to be particularly effective at dis- 
tinguishing between a radial or meridional field from an az- 
imuthal field. They also found crosstalk from radial to merid- 
ional field components and v ice versa, especially at low lati- 
tudes. Donati & Brownl (1 19971) used an even simpler weak-field 
Gaussian analytical treatment of the line profiles, which does not 
allow realistic computation of the polarization prof iles for mag- 
netic s tars with temperature inhom ogeneities. Both lBrown et al.l 
dl991l) and [Donati & Brown) d 19971) constrained their inversions 
by the Maximum Entropy method. 

Due to weakness of surface magnetic fields on cool stars, 
they are usually observed in Stokes / and V parameters only. 
Using a n independent ZDI code empl oying Tikhonov regular- 
ization, Kochukhov & Piskunov (2002) investigated the differ- 
ence in using all four Stokes parameters for magnetic inversions 
compared to using only circular polarization. In agreement with 
Donati & Brown (1997), they found a strong crosstalk between 
radial and meridional field components at low latitudes, which 
can be alleviated if inversions are base d on the full Stokes vec- 
tor ob servations. However, the tests by [Kochukhov & Piskunov! 
(2002) targeted early-type Ap stars with chemical spots and did 
not incorporate temperature inhomogeneities in magnetic field 
inversions nor did they try to simultaneously reconstruct the 
magnetic field and temperature maps. 

Summarizing, the conventional method of mapping mag- 
netic fields in cool stars is to use Stokes V line profiles only, 
ignoring local temperature variations which can be derived from 
simultaneous analysis of Stokes /. Previous ZDI numerical ex- 
periments also mainly focused on purely magnetic or purely 
temperature spots. Yet, from general theoretical considerations 
and observations of the Sun we expect a strong correlation of 
the temperature and magnetic features. It is not clear just by 
how much this correlation can undermine magnetic maps recon- 
structed from Stokes V alone. 

The goal of our paper is to perform a series of numerical 
tests to verify validity of ZDI maps obtained for stars with tem- 
perature spots and to investigate if the quality of magnetic map- 
ping would increase if temperature inhomogeneities are taken 
into account in a self-consistent and physically realistic manner. 
The inversion code we have used is suitable for this task since 
it computes the Stokes line profiles self-consistently, taking into 
account both the local Zeeman effect and the local temperature. 

Our paper is structured as follows. We describe our ZDI code 
in Sect. [2] and discuss the setup for numerical experiments in 
Sect. [3] We then present inversions from the simulated data in 
Sect. HI evaluate these results in Sect. |5j and summarize main 
conclusions of our investigation in Sect. [6] 

2. Zeeman Doppler Imaging code 

In this paper we carried out ZDI numerical experiments us- 
ing our new inversion code Invers13. This magnetic imag- 
ing software has been developed from the InversIO cod e 
(Pisk unov & Kochukhovl [20021: IKochukhov & Piskunov! 12002b . 
previously applied for the reconstruction of magnetic field 
topologies and distributions of chemical spots on Ap stars (e.g., 
IKochukhov & Waddl2010t Oiftinger e t~al1l2010l) . Invers13 in- 
corporates the full numerical treatment of the polarized radia- 
tive transfer in a realistic stellar model atmosphere. This allows 
us to calculate the spectra in all four Stokes parameters simul- 
taneously and self-consistently, without relying on any of the 
common simplifying assumptions, such as Gaussian profiles in 



the weak-field limit or the Milne-Eddington atmosphere. Of par- 
ticular importance for the magnetic mapping of cool stars is a 
self-consistent treatment of the temperature and magnetic stellar 
surface structures: the temperature variations are accounted for 
in computing the polarization profiles and, at the same time, the 
Zeeman splitting of spectral lines is incorporated in theoretical 
Stokes / spectra. The approach implemented in Invers13 is su- 
perior to most other applications of ZDI, which systematically 
neglect the effects of temperature spots in the magnetic inver- 
sions for cool active stars. 

Calculation of the local Stokes parameter profiles by 
Invers13 is b ased on the spectra l line lists extracted from the 
vald database (iKupka et al.l ll999). The model atmospheres for a 
given metallicity, surface gravity and a r ange of effective temper - 
atures are adopted from the marcs grid (Gustafs son et al.l l2008). 
Invers13 calculates the Stokes IQUV profiles and continuum in- 
tensities for an arbitrary effective temperature of a given stellar 
surface element by interpolating between the sets of model spec- 
tra corresponding to the three nearest points in the model atmo- 
sphere grid. This allows an accurate semi-analytical computa- 
tion of the derivatives with respect to temperature. The deriva- 
tives with respect to the three magnetic field vector components 
are evaluated numerically, using a simple one-sided difference 
scheme. 

The local spectra are convolved with a Gaussian profile 
to take into account instrumental and radial-tangential macro- 
turbulent broadening. Then, the Stokes spectra are Doppler- 
shifted, interpolated on the wavelength grid of the observed 
spectra and summed for each rotational phase. The resulting 
disk-integrated Stokes parameter profiles are normalized by the 
phase-dependent continuum flux and compared with observa- 
tions. In calculating the goodness of fit it is essential to balance 
the contributions of each Stokes parameter by scaling the respec- 
tive chi-square terms by the inverse of the mean amplitude of the 
corresponding Stokes parameter. This usually implies a 10 to 30 
times higher weight for the Stokes V profiles compared to Stokes 

The ZDI with only Stokes / and V spectra available for 
cool active stars is an intrinsically ill-posed problem, requir- 
ing the u se of regularizatio n to reach a stable and uniqu e 
solution (Bro wn et af] Il99lb IPiskunov & Kochukhovl 120021) . 
Reconstruction of temperature in Invers13 is stabilized using 
the Tikhonov method, similar to many prev ious DI studies of 
cool active stars (e.g., IPiskunov & Ricdll993h . Tikhonov regu- 
larization ensures that the code converges on a surface distri- 
bution with a minimum contrast between adjacent surface ele- 
ments. For the magnetic field reconstruction, Invers13 provides 
a choice between either directly mapping the three magnetic vec- 
tor components and applying the Tikhonov regularization indi- 
vidually to rad ial, meridional and azimutha l magnetic maps as 
described by Pisk unov & Kochukhovl d2002l) or expanding the 
field into a spherica l harmonic series following the approach by 
IDonati et al.l J200 6!). We used the former regularization method 
in the ZDI numerical experiments presented in this paper be- 
cause our goal was to investigate reliability of reconstruction of 
isolated small-scale magnetic features. 

The iterative adjustment of the surface magnetic and tem- 
perature maps to match available observations and satisfy regu- 
larization constraints is accomplished b y means of the modified 
Leven berg-Marquardt optimization (Pis kunov & Kochukhovl 
2002). This algorithm enables convergence in typically 10-20 
iterations starting from a homogeneous temperature distri- 
bution and zero magnetic field. Invers13 code is optimized 
for execution on parallel computers using MPI libraries. For 
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Table 1. Configuration of the magnetic and temperature spots adopted for ZDI numerical experiments. 



Latitude 


Longitude 


Spot radius (R*) 


B r (kG) 


B a (kG) 


B, (kG) 


T s (kK) 








0.40 


1.0 








4.75 


20 


90 


0.40 


1.0 








4.75 


40 


180 


0.40 


1.0 








4.75 


60 


270 


0.40 


1.0 








4.75 








0.25 


2.0 








3.75 


20 


90 


0.25 


2.0 








3.75 


40 


180 


0.25 


2.0 








3.75 


60 


270 


0.25 


2.0 








3.75 



Notes. A similar setup was used for the experiments with meridional and azimuthal fields inside cool spots. 



this particular study, we were running the code on a 16-CPU 
HP V-class 2500 server. Further details on the numerical and 
computational methods employed by Invers13 can be found in 
IKochukhov et all d2012l) . 

3. Setup for numerical experiments 

For the calculations in this paper we used a grid consisting of 
1176 surface elements with roughly equal areas. Each element 
was assigned a specific value of magnetic field strength and tem- 
perature. The inclination angle i of the star was set to 50°, the 
projected rotational velocity was v e sin/ = 40 kms~' and the 
photospheric temperature was 5750 K. We created four circu- 
lar magnetic spots at different latitudes and longitudes. By do- 
ing so we could investigate the quality of mapping as a func- 
tion of spot position, similar to the experiments in previous stud- 
ies dDonati &Brownlll997l: IKochukhov & Piskunovll2002l) . All 
spots were identical in shape and field strength, with a differ- 
ence of 1 kG between the umbra and penumbra. The spots were 
circular, with a radius of 0.25/?+ for the umbra and 0.4/?+ for the 
penumbra. All magnetic spots also had the same field orienta- 
tion: either radial, meridional or azimuthal. The surface outside 
the spots was assumed to be non-magnetic. 

We used two different temperature distributions. The first 
was homogeneous, i.e. the entire surface including the area cov- 
ered by magnetic spots had a temperature equal to the photo- 
spheric temperature. This configuration of homogeneous tem- 
perature and magnetic spots, used in the first set of experi- 
ments, represents a reference for the subsequent reconstructions 
of fields inside cool spots. It may also correspond to a real situ- 
ation for active stars with a hypothetical non-solar dynamo pro- 
ducing temperature structures not associated with particularly 
strong magnetic fields, or for low-activity stars for which tem- 
perature inhomogeneities are u nresolved and on ly large-scale 
magnetic fields are detected (e.p. lPetit et al.ll2008l) . 

The second temperature distribution consisted of four cool 
temperature spots with the exact same shapes and positions as 
the magnetic spots. The umbra and penumbra was, respectively, 
2000 K and 1000 K cooler than the surrounding photosphere. 
This configuration of temperature and magnetic spots was used 
in the second and third sets of experiments. A complete descrip- 
tion of the surface structure setup for the calculation using radial 
magnetic spots can be found in Table [TJ 

For all th e experiments we used a grid of 22 marcs model 
atmospheres dGustafsson et al.l 120081) with the surface gravity 
log g = 4.5 and microturbulence = 2 km s _1 . Since the temper- 
ature of the umbra was set to 3750 K and the photospheric tem- 
perature was 5750 K, the marcs model atmosphere grid had a 
range in T e ff between 3000 K and 6750 K in order to cover these 
temperatures. We used three magnetically sensitive Fei lines at 



5497.5, 5501.5 and 5506.8 A. To make our calculations realis- 
tic, we also included seven metal lines blending the Fe i lines of 
interest. We extract ed atomic line data from the vald data base 
( Kup ka et al.l 1 1 9991) to get information necessary for spectrum 
synthesis and to obtain Lande factors and quantum numbers re- 
quired for Zeeman splitting calculations. These data were then 
used as input for forward calculations in order to obtain the ro- 
tationally modulated line profiles of either Stokes / and V, or all 
four Stokes parameters, at 20 equally spaced rotational phases. 
Profiles were convolved with a Gaussian function to simulate in- 
strumental resolution of A/AA = 65000, but no noise was added 
to the profiles since we wanted to test the performance of ZDI 
under ideal conditions. That is also why we chose to use such a 
good phase coverage and optimal values for v e sin i and i. Other 
ZDI studies have investigated the effect of different phase cov- 
erage, chJJerent^_sn^ di^ noise 
level (iDonati & Browdll997tlkochukhov & Piskunovll2002l) . 

The resulting line profiles were then used as an input for in- 
verse calculations to recover the surface maps. In experiment 1 
and 2 we first used only Stokes V to reconstruct the magnetic 
field maps for each of the three magnetic field vector compo- 
nents. Then we used both Stokes / and V to reconstruct the mag- 
netic field and temperature map. In the last experiment we used 
Stokes /, Q, U and V parameters for magnetic and temperature 
mapping. Below we present results of the inversions for these 
three sets of tests. 



4. Results 

In the following sections we discuss results of the simulated ZDI 
inversions for purely magnetic spots using Stokes IV, for mag- 
netic field concentrations coinciding with low-temperature spots 
by first using Stokes IV, and, finally, by using all four Stokes pa- 
rameters for the inversions. In each case we reconstructed mag- 
netic and temperature maps for different magnetic field vector 
orientations. We assess these inversion results by the visual com- 
parison of the true and reconstructed surface distributions rep- 
resented in rectangular projection maps. We also evaluate the 
inversions numerically, by computing by how much the mag- 
netic field strength is underestimated in the central regions of 
magnetic spots at different latitudes. Since the Tikhonov regu- 
larization applied by our code always smoothes the horizontal 
gradients, leading to some underestimation of both temperature 
and magnetic field contrast for small spots, we also calculate a 
complementary measure of the total magnetic field modulus as 
well as the total radial, meridional, azimuthal fields by integrat- 
ing corresponding distributions over the entire stellar surface. 
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Fig. 1. Initial and reconstructed magnetic field maps for the inversions with a fixed homogeneous temperature. The upper panel in 
each pair of rows corresponds to the true maps while the lower panel shows the reconstructed maps. These inversions were carried 
out using only Stokes V and assuming a constant temperature. 



Table 2. Underestimation of field strength for the inversion with 
a homogeneous temperature. 



Used Stokes Field Spot 1 Spot 2 Spot 3 Spot 4 

parameters component lat. 60° lat. 40° lat. 20° lat. 0° 

V ~B r 21% 31% 49% 78%" 
Btot 21% 31% 40% 71% 

V B m 42% 69% 58% 70% 
gto, 42% 68% 52% 64% 

V B a 30% 40% 28% 56% 
gtoi 30% 40% 28% 55% 

IV fi r 10% 12% 12% 95% 
9% 12% 11% 22% 

IV B m 6% 12% 6% 20% 
But 6% 11% 6% 20% 

IV S a 7% 8% 7% 13% 

B tot 7% 8% 7% 13% 



Notes. The compared surface elements are taken from the center of each 
spot. 



4.7. Homogeneous temperature distribution 

The maps resulting from the inversions using line profiles simu- 
lated for the homogeneous temperature distribution are shown 
in Figs. Q] and [2] There are clear differences between the re- 



constructed magnetic field maps in Fig. Q] and those in Fig. [2] 
even though the only difference between these inversions was 
that the code was allowed to vary temperature in the latter case. 
The most prominent difference between the two sets of maps is 
the field strength. A closer inspection was made by comparing 
the field values of specific surface elements of the star in the re- 
constructed maps to the same elements in the true maps. Four 
surface elements from the center of each spot were chosen and 
an average discrepancy value was calculated. These results are 
reported in Table[2]both for the field modulus and for the specific 
field component (a large difference between the two values indi- 
cates significant crosstalk). The spot positions are however quite 
accurately reconstructed in both cases even though some spots, 
especially in Fig.Q] appear a bit elongated. However, the center 
of the spots, in terms of the highest magnetic field strength, does 
not suffer from any shifts in longitude nor in latitude. 

The field strengths of individual surface zones are only un- 
derestimated by less than 10% when both Stokes / and V were 
used for the inversion. When temperature is kept constant the 
field strength underestimation of individual components varies 
between 21% and 78%. We also integrated each field compo- 
nent and the total field strength over the entire stellar surface. 
The resulting values can be found in Table [3] 

By inspecting the reconstructed temperature maps shown in 
Fig.|2]it appears that the areas covered by magnetic spots contain 
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Fig. 2. Same as for Fig.Q]but here both Stokes / and V were used and possible temperature inhomogeneities were taken into account 
in self-consistent ZDI inversions. 



Table 3. Integrated field values for the inversions with a homogeneous temperature. 



Used Stokes 


Field 


True 


Reconstructed 


True 


Reconstructed 


True 


Reconstructed 


parameters 


component 


map (kGrad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


V 


B r 


2.945 


2.682 





0.340 





0.047 




B m 





0.798 


2.945 


2.997 





0.030 




B, 





-0.084 





-0.052 


2.945 


3.616 




B to t 


2.945 


3.574 


2.945 


3.550 


2.945 


4.117 


IV 


B T 


2.945 


3.828 





-0.687 





0.265 




B m 





0.497 


2.945 


4.065 





-0.305 




B, 





-0.099 





-0.009 


2.945 


3.408 




B tot 


2.945 


5.360 


2.945 


4.524 


2.945 


4.378 



some temperature variations. However, the full range of temper- 
ature variation is only about 50-100 K. The temperature spots 
appear slightly hotter than the rest of the stellar surface. 

When temperature is assumed constant the reconstructed ra- 
dial field spots (see Fig. [TJ are all underestimated in strength. 
According to Table [2] the underestimation increases with de- 
creasing latitude. This is also true for the radial spots in Fig. [2] 
although the underestimation is not as severe except for the low- 
est latitude spot which is instead interpreted as a meridional spot 
by the inversion code. The total field strength of that spot is how- 
ever only underestimated by 22%, as can be seen in Table[2] The 
total field strength and hence also the radial field strength inte- 
grated over the entire stellar surface is 2.945 kGrad 2 in the true 
map, see Table [3] When both magnetic field and temperature are 



reconstructed, the total field strength and radial field strength are 
higher, 5.360 kGrad 2 and 3.828 kGrad 2 respectively. This can 
also be seen in Fig. [2] since almost the entire map is covered 
by a radial field including three spots which are only underesti- 
mated with about 10%. When temperature is kept constant, the 
total field strength is slightly overestimated while the radial field 
strength instead is slightly underestimated. The meridional field 
value, which should be zero, is instead almost 800 G rad 2 , and 
this can also be seen as crosstalk between radial and meridional 
fields in Fig. Q] especially at low latitudes. The two lowest lati- 
tude spots has an almost as strong meridional component as the 
radial one. The increasing crosstalk with decreasing latitude is 
also evident from the values in Table [2] The difference in field 
underestimation between the radial field component and the to- 
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Fig. 3. Same as for Fig.[T]but for the inversions with the magnetic spots coinciding with low-temperature spots. These inversions 
were carried out using only Stokes V and assuming a constant temperature. 



tal field component systematically increases with decreasing lat- 
itude. 

The biggest discrepancy between the true and recovered field 
maps in Figs. Q] and [2] can t> e found for meridional field spots. 
By comparing the strength of the same spot for each of the three 
field directions in Table |2]it is evident that the meridional spots 
are the weakest for three out of the four spots when only Stokes 
V is used for inversions. The underestimation of the meridional 
field component varies between 49% and 70% in this case. The 
situation is almost the opposite when both Stokes / and V are 
used for inversions since now two of the meridional spots are 
the strongest since they are only underestimated by 6-20%. The 
integrated total and meridional field values are however overes- 
timated in both reconstructions. Crosstalk can also be seen and 
once again it is strongest between meridional and radial compo- 
nents at low latitudes. The overall structure of the spots in the 
meridional maps in Figs. Q] and [2] are somewhat similar in the 
sense that the two high latitude spots are connected, and the two 
low latitude spots are spread out towards lower latitudes. Even 
so, the spots in Fig. [2] are better defined since there is a larger 
difference in field strength between the actual spots and the field 
outside the spots. 

The reconstructed azimuthal magnetic field spots suffer least 
from crosstalk but are however quite similar to the meridional 
field maps in terms of structure. The two high-latitude spots are 
overlapping while the two low-latitude spots spread out towards 



lower latitudes. The spots in Fig.|2]are still more clearly defined 
than for the corresponding inversion using Stokes V alone. Once 
again the integrated field strengths of the reconstructed maps ex- 
ceed those of the true map. The field strength is somewhere in 
between the values for the radial and meridional spots at high 
latitudes, but higher at low latitudes. 

Table 4. Underestimation of field strength for the inversions with 
an inhomogeneous temperature distribution. 



Used Stokes Field Spot 1 Spot 2 Spot 3 Spot 4 

parameters component lat. 60° lat. 40° lat. 20° lat. 0° 

V B~ r 82% 85% 89% 94%" 
B tot 82% 85% 87% 92% 

V B m 87% 95% 94% 94% 
g, ot 86% 94% 90% 91% 

V B., 83% 85% 85% 91% 
ff to , 83% 85% 85% 91% 

IV Br 46% 45% 60% 76% 
B tot 45% 45% 55% 67% 

IV B m 70% 84% 70% 73% 
B tot 69% 81% 64% 67% 

IV B., 46% 53% 47% 60% 

Btot 46% 53% 47% 60% 



Notes. The compared surface elements are taken from the center of each 
spot. 
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Fig. 4. Same as for Fig. [3] but here both Stokes / and V were used and temperature inhomogeneities were taken into account in 
self-consistent ZDI inversions. 



4.2. Inhomogeneous temperature distribution 

The results of the inversions for magnetic features coinciding 
with low-temperature spots are presented in Figs.[3]and|4] As can 
be seen from the reconsttucted maps shown in Fig. [3] the mag- 
netic field strength is strongly underestimated for all three field 
components when temperature variations are ignored. By com- 
paring the field values of specific surface elements in the same 
way as in the first experiment, we inferred that the field strengths 
are underestimated by as much as 80-95% (see Table|4j>. On the 
other hand, the temperature reconstruction in the inversion using 
both Stokes / and V once again proved to be very accurate. There 
is a tendency for a slight overestimation, but only by approxi- 
mately 250 K. As in the first experiment, the field strengths of 
the reconstructed radial spots shown in Fig. [3] seem to decrease 
with decreasing latitude. When temperature variations are taken 
into account (Fig. |4| the two high-latitude spots have similar 
sttengths while the two low-latitude spots are distorted similarly 
to the first experiment. Once again there is crosstalk primarily 
between the radial and meridional field components at low lati- 
tudes in both cases. This can be seen directly in the reconstructed 
maps, but also in Table [4] since once again the underestimation 
of the total field strength is less compared to the underestima- 
tion of the radial field strength at low latitudes. Otherwise there 
are not many similarities between the two radial field maps. The 
underestimation of the radial field strengths when only Stokes 
V is used for inversion varies between 82% and 94%. The un- 



derestimation is not as severe when temperature is taken into ac- 
count and varies between 45% and 76% instead. The integrated 
field values, which can be found in Table [5] show that the radial 
field component and the total field strength are almost twice as 
large when Stokes / and V are used for inversions compared to 
when only Stokes V is used. Despite this, the meridional field 
sttengths, which should be zero when the true spots are purely 
radial, are comparable in strength, and so are the azimuthal field 
strengths. In other words, even though the radial field strength is 
so much weaker when temperature variations are ignored, there 
is still almost as much crosstalk as when temperature is taken 
into account. The spot geometry is quite accurately reproduced 
in Fig. [4] (self-consistent magnetic and temperature inversion 
from both Stokes / and V) while it is almost impossible to relate 
to the original maps in Fig. [3] (magnetic mapping from Stokes V 
alone). 

Comparing the quality of all the reconstructions in all the ex- 
periments, the worst one can be found in the experiment for the 
meridional spots when temperature variations are not taken into 
account. The meridional field strength is underestimated by as 
much as 87% for one spot and by 94-95% for the other three 
spots. The situation is somewhat improved when temperature 
is taken into account since the field is underestimated by 70- 
84% instead. The integrated total and meridional field strengths 
are almost three times higher when inversions are made using 
both Stokes / and V compared to when only Stokes V is used. 
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Table 5. Integrated field values for the inversions with an inhomogeneous temperature distribution. 



Used Stokes 


Field 


True 


Reconstructed 


True 


Reconstructed 


True 


Reconstructed 


parameters 


component 


map (kG rad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


V 




2.945 


1.243 





0.235 





0.072 




B m 





0.423 


2.945 


0.553 





0.003 




B a 





-0.090 





-0.016 


2.945 


1.313 




Slot 


2.945 


1.479 


2.945 


0.864 


2.945 


1.391 


IV 


Br 


2.945 


2.159 





0.647 





-0.091 




B m 





0.474 


2.945 


1.645 





0.184 




B a 





-0.078 





-0.214 


2.945 


2.233 




Bioi 


2.945 


2.994 


2.945 


2.539 


2.945 


2.721 




Radial field 


Meridional field 


Azimuthal field 


Temperature 





(i — 

o.o — 

-1.0 



(i — 

0.0 

-1.0 

(>-■ 



90 ~ 
60 - 
•3 30 : 

f o- 

3 -30 - 

B[kO] 

90 r~ 

60 - 
% 30 - 
I 0: 
3 -30 - 
-60 " " 



-90 L 

























f W * • 











T[kK] 



60 120 180 240 300 60 120 180 240 300 60 120 180 240 300 60 120 180 240 300 360 
Longitude Longitude Longitude Longitude 



Radial field 



Meridional field 



Azimuthal field 



Temperature 



T[kK] 



60 120 180 240 300 
Longitude 



60 120 180 240 300 60 120 180 240 300 60 120 180 240 300 360 
Longitude Longitude Longitude 



Meridional field 



Azimuthal field 



Temperature 




T[kK] 



60 120 180 240 300 60 120 180 240 300 60 120 180 240 300 60 120 180 240 300 360 
Longitude Longitude Longitude Longitude 



Fig. 5. Same as for Fig. |4]but here all four Stokes parameters were used in self-consistent ZDI inversions. 



Table 6. Underestimation of field strength for the self-consistent 
inversions using all four Stokes parameters. 



Used Stokes 


Field 


Spot 1 


Spot 2 


Spot 3 


Spot 4 


parameters 


component 


lat. 60° 


lat. 40° 


lat. 20° 


lat. 0° 


IQUV 


B r 


33% 


36% 


41% 


49% 




B to t 


32% 


36% 


40% 


49% 


IQUV 


B m 


31% 


44% 


41% 


62% 




Biot 


31% 


44% 


41% 


61% 


IQUV 


B a 


34% 


37% 


36% 


48% 




Btot 


34% 


37% 


36% 


43% 



Notes. The compared surface elements are taken from the center of each 
spot. 
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Table 7. Integrated field values when all four Stokes parameters are used for inversions. 



Used Stokes 


Field 


True 


Reconstructed 


True 


Reconstructed 


True 


Reconstructed 


parameters 


component 


map (kG rad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


map (kGrad 2 ) 


IQUV 


B t 


2.945 


2.480 





0.099 





-0.017 




B m 





0.054 


2.945 


2.363 





-0.096 




Ba 





-0.001 





0.019 


2.945 


3.029 




Slot 


2.945 


3.157 


2.945 


2.802 


2.945 


3.458 



However, there is still a clear crosstalk between the meridional 
and radial field components (see Fig. |4j>. 

The reconstruction of the azimuthal field spots follows the 
same pattern as for the other two field directions in the sense 
that recovery of the spot strength and geometry are significantly 
improved when temperature variations are taken into account. 
The azimuthal field strengths are underestimated by 46-60% 
compared to 83-91% for the inversion neglecting temperature 
spots. The integrated total azimuthal field strength is twice as 
large. Similar to the first experiment, the azimuthal spots are 
most robust and suffer very little from crosstalk. Furthermore, 
the strengths of the two lowest latitude spots are higher com- 
pared to the corresponding radial and meridional spots. 

4.3. All four Stokes parameters 

The maps obtained in the inversions using all four Stokes param- 
eters are shown in Fig. [5] The field strength is underestimated by 
31% to 62% (see Table |6j, but the spot geometry is very ac- 
curately reconstructed, with the only exception being the low- 
est latitude azimuthal spot. Almost none of the crosstalk seen in 
the first two experiments is visible in the full Stokes vector in- 
version. The temperature is once again correctly reproduced to 
within 200 K in general. 

As in the other two experiments, the radial field strength 
underestimation is systematically increasing from 33% for the 
highest-latitude spot to 49% for the lowest-latitude spot. The to- 
tal field strength is higher than the true value, but the radial field 
strength is lower, see Table [7] The spot geometry and positions 
are accurately recovered and the map shows little spurious mag- 
netic flux outside the spots. No significant crosstalk is detected 
and it appears that the ability of the inversion code to distinguish 
the radial and meridional component is now as good as it was 
for the radial and azimuthal components in the inversions with 
circular polarization alone. 

The geometry of meridional spots is correctly recovered and 
their strengths are underestimated by 31-62%. Here the inte- 
grated total field and meridional field values are slightly under- 
estimated. 

The azimuthal spots have similar strengths to the radial spots 
since they are underestimated by 34^8%. The lowest-latitude 
spot does however seem to spread over lower latitudes making 
the reconstructed spot topology less well-defined compared to 
the other two maps. This can also be seen in the integrated field 
values, which are slightly overestimated even though the field 
strengths of the centers of spots are underestimated. Once again, 
the reconstructed field strengths inside the two low-latitude spots 
are higher for the azimuthal field spots compared to the same 
experiments with radial and meridional field spots. 

5. Discussion 

In this study we carried out numerical ZDI simulations by in- 
verting spectra generated using a standard grid of stellar model 



atmospheres. In other words, following numerous previous con- 
ventional temperature DI studies, we have assumed that the 
thermal structure of a starspot can be approximated by a sim- 
ple radiative equilibrium model atmosphere. Studies of the Sun 
demonstrated validit y of such an approximation at least fo r the 
sunspot interiors (see Solanki & Unruh 1998; Solanki 2003, and 
references therein). Until realistic magnetic starspot models be- 
come widely available, this is the only feasible approach for 
the model-atmosphere-based active-star DI. In any case, we are 
mostly concerned with including in the magnetic mapping a 
dominant effect of the intensity contrast due to presence of hot 
and cool surface features. The temperature dependence of the 
local Stokes / profiles is the second most important factor, prop- 
erly incorporated in our calculations. On top of that, possible de- 
viations of the T - t relation from the predictions of hydrostatic 
model atmospheres might induce line shape changes of Stokes 
parameters. But, even if present, this effect is going to be by far 
less important than the scaling of the polarization profiles by the 
local continuum intensity and by the factor reflecting tempera- 
ture sensitivity of specific spectral lines. 

One of the unexpected outcomes of our numerical experi- 
ments is the difference of magnetic maps inferred from Stokes 
IV and from V alone in the absence of temperature spots. When 
the temperature is homogeneous one might expect that the re- 
sulting magnetic field maps should not depend on whether tem- 
perature was reconstructed since in both cases the code ends with 
a nearly identical temperature distribution lacking significant in- 
homogeneities. But the magnetic maps presented in Figs.[TJand[2] 
are not identical although the temperature is the same to within 
50-100 K. The explanation to this discrepancy might be that in- 
formation about the magnetic field is also contained in Stokes / 
which is made available and taken into account when Stokes / is 
used for simultaneous temperature mapping. In order to test this 
further, we performed another inverse calculation of the stellar 
surface with magnetic spots and a homogeneous temperature us- 
ing both Stokes / and V, but excluding the Zeeman effect in the 
forward calculations of Stokes /. The magnetic maps obtained 
in this case resemble more closely the magnetic distributions 
shown in Fig. Q~| which were reconstructed using only Stokes 
V, rather than the ones illustrated in Fig. [2] 

We found that for our assumed magnetic spot configuration 
Stokes / is indeed affected by the magnetic field. The Zeeman ef- 
fect leads to splitting and broadening of local spectral lines and 
this shows up as a broadening and deepening of the correspond- 
ing disk-integrated Stokes / profile. In Fig. [6] we have plotted 
all Stokes / profiles from the first experiment, i.e. when the stel- 
lar surface had magnetic spots but a homogeneous temperature 
together with another set of Stokes / profiles of a stellar sur- 
face with the same homogeneous temperature but without any 
magnetic spots. There is a clear difference between the two sets 
of profiles since the one where magnetic spots were included is 
rotationally modulated. Thus, not only does temperature varia- 
tions have to be taken into account when magnetic mapping is 
done, but a magnetic field also have to be taken into account 
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Fig. 6. Stokes / profiles from a stellar surface with homogeneous 
temperature. The black solid lines represents the Stokes / pro- 
files when the stellar surface contained magnetic spots and the 
red dashed lines represents Stokes / without any magnetic spots. 

when doing temperature mapping. That is certainly one of the 
reasons to why the temperature reconstruction is accurate in our 
experiments. If only temperature is reconstructed using Stokes 
/ in the case with magnetic spots on a surface with homoge- 
neous temperature, the spurious temperature variations reach to 
about 500 K locally. This should be compared to a difference of 
100 K when Stokes V is included for simultaneous temperature 
and magnetic field reconstruction as in the first experiment. 

We also performed the same test for the surface distribution 
used in the second experiment, i.e. the map with cool magnetic 
spots. When only Stokes / was used to reconstruct the temper- 
ature map ignoring magnetic field, the temperature variations 
were found to be close to the temperature distribution recon- 
structed self-consistently using both Stokes / and V. This is, per- 
haps, not so unexpected since the influence of the Zeeman effect 
on Stokes / is significantly smaller than the distortions caused 
by the cool temperature spots. These results imply that the con- 
ventional temperature DI, which does not take the Zeeman ef- 
fect into account, is not necessarily incorrect, at least for the 
magnetic field strengths comparable to those investigated in our 
study. 
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Fig. 7. Line profiles from forward and inverse calculations. The 
black and red lines correspond to the line profiles from forward 
and inverse calculations respectively of a stellar surface with 
magnetic spots and cool temperature spots. The blue lines repre- 
sents line profiles from forward calculations of a stellar surface 
with magnetic spots and a homogeneous temperature distribu- 
tion. 
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Perhaps it is not so surprising that the magnetic field recon- 
struction is so poor when temperature variations are ignored. The 
amplitude of the Stokes V line profiles depends on both temper- 
ature and magnetic field strength. This has also been discussed 
before, for example by Don ati & Brow r3( fl99l . Instead of si- 
multaneously reconstructing temperature, they are using a field 
strength of 500 G which they say corresponds to a 2 kG field, in 
Stokes V, with 50% filling factors and 50% temperature-induced 
signal dilution. When temperature is not reconstructed simul- 
taneously with magnetic field, the temperature of the magnetic 
spots is assumed to be the same as the photospheric temperature. 
The low amplitude of the polarized line profile is then incor- 
rectly interpreted as if the magnetic field strength is low. When 
temperature is taken into account, as in our magnetic inversion 
technique, the code adjusts the field strength according to the ap- 
propriate local temperature. The low amplitude of the polarized 
profile is therefore not misinterpreted in terms of low magnetic 
field strength. To illustrate the striking difference in amplitude 
of magnetic signal associated with cool spots compared to mag- 
netic spots at photospheric temperature, we show the respective 
Stokes V profiles in Fig. [7] 

Furthermore, besides a general major underestimate of the 
magnetic flux in the ZDI without accounting for temperature 
spots, our experiments suggest that magnetic field topology also 
cannot be correctly reconstructed with such an inconsistent ap- 
proach. Comparison of Figs. [3] and [4] clearly shows that mag- 
netic reconstruction ignoring temperature spots does not yield a 
faint copy of the correct map. Instead, we get highly distorted 
magnetic configurations, which in many cases bear little resem- 
blance to the input magnetic maps. Thus, contrary to the assump- 
tion of the ZD I studies which treat Stokes / and V fully inde- 
pendently (e.g. lDonati & Brownlll997uDonatill999l) . the lack of 
self-consistency between interpretation of Stokes / and V actu- 
ally jeopardizes all aspects of magnetic mapping. 

We have verified that these inversion results are not due to a 
large difference in the quality of profile fits for different inver- 
sions. Fig. [8] shows the typical fits to line profiles correspond- 
ing to different experiments. Evidently, the difference is not 
that large. The worst fits are, not surprisingly, found for fixed- 
temperature fitting of the Stokes V profiles associated with an 
inhomogeneous temperature distribution. In this case, a single- 
component atmosphere assumed by the code is physically incor- 
rect and cannot reproduce all details of the circular polarization 
profile by adjusting the field strength and orientation. But even 
in this case the fits are not grossly incorrect in the sense that 
the profile amplitudes are not severely underestimated. What can 
also be seen from Fig.[8]is that the quality of the fit does not sig- 
nificantly depend on the field direction. 

We found that even in the simultaneous inversion using both 
Stokes / and V magnetic reconstruction suffers a number of 
crosstalks. In particular, the crosstalk between the radial and 
meridional field components is prominent and can be explained 
by the similarity in the corresponding polarized line profiles 
dDonati & Brownlll997l) . Circularly polarized light is diagnos- 
ing the line-of-sight component of the magnetic field vector. A 
magnetic spot with a meridional field will therefore produce an 
anti-symmetric Stokes V signature, quite similar to that of a ra- 
dial field spot, since the projection of the magnetic field vector 
on the line-of-sight will not change sign with rotation phase for 
a given spot. An azimuthal field will, however, produce a differ- 
ent dynamical Stokes V signature since the projected magnetic 
field vector will change sign in the center of the stellar disk as 
the spot moves relative to the observer due to stellar rotation. In 
Fig [8]Stokes V line profiles taken from the same rotational phase 



Radial field 



Meridional field 



Azimuthal field 






(a) Line profiles when temperature is homogeneous and only Stokes 
V is used for inversions. 






(b) Line profiles when temperature is homogeneous and both Stokes 
/ and V are used for inversions. 





(c) Line profiles when temperature is inhomogeneous and only 
Stokes V is used for inversions. 






(d) Line profiles when temperature is inhomogeneous and both 
Stokes / and V are used for inversions. 
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(e) Line profiles when temperature is inhomogeneous and all four 
Stokes parameters are used for inversions. 

Fig. 8. Example of Stokes V profile fits for the rotational phase 
0.1 for each of the inversion tests with radial field. Black lines 
represent the true line profiles produced with forward calculation 
and the red lines show corresponding fits by Invers13 code. 



can be seen. The profiles corresponding to radial and meridional 
field spots are similar in shape while the profiles representing 
azimuthal field spots are very different compared to the other 
two. 

The crosstalk between radial and meridional fields is 
strongest at low latitudes. This feature is easy to explain because 
in this case magnetic spots will be observed from a limited range 
of viewing angles. In general, the magnitude of the line-of-sight 
projection of the magnetic field vector for a spot with a constant 
radial field will change with phase while the magnitude of a pro- 
jected field vector for a spot with a constant meridional field will 
be almost constant with rotational phase. This allows the code 
to distinguish radial and meridional fields at higher latitudes 
which are observable over a larger range of phases. However, 
low-latitude spots are seen from approximately the same aspect 
angle, resulting in approximately constant projected field vector. 
This makes it difficult to distinguish radial and meridional fields. 

The meridional field is the most difficult field orientation to 
reproduce. On the one hand, as already mentioned, a meridional 
field might be interpreted as a radial field because of the similar- 
ity in the respective Stokes V profile shapes. On the other hand, 
a field vector of a radial or an azimuthal field will at some ro- 
tational phase point more or less straight towards the observer, 
depending on the stellar inclination and phase of course, mak- 
ing the projected field vector almost as large as the local field 
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vector. A meridional field vector will, however, always be al- 
most perpendicular to the line-of-sight, making the line-of-sight 
magnetic vector projection small and hence resulting in a small 
amplitude of the circular polarization profile. These aspects con- 
spire to make recovery of the meridional field very challenging. 

When all four Stokes parameters are used, most of the 
crosstalk is removed. By including linear polarization, which is 
sensitive to the transverse component of the magnetic field vec- 
tor, a lot more information about the orientation of the field can 
be retrieved, as alread y demonstrated by the ZDI ex periments 
for magnetic Ap stars dKochukhov & Piskunovl|200 2). By com- 
paring Fig.[3]to Fig. [5] which have the same spot set-up, one can 
see that the reconstructed fields are generally stronger when all 
four Stokes parameters are used, especially for the meridional 
field. Even though the line-of-sight component of a meridional 
field vector is small, the perpendicular component is large and 
this is detected by the linear polarization. 



6. Conclusions 

We have carried out a series of numerical experiments with 
Zeeman Doppler imaging, aiming to explore reliability of re- 
construction of different magnetic field geometries and testing to 
what extent independent reconstruction of temperature and mag- 
netic spots is meaningful for stars with high-contrast temperature 
inhomogeneities. We performed these tests under ideal condi- 
tions, adopting optimal stellar parameters in terms of v e sin i and 
i, using an excellent phase coverage and adding no noise to the 
simulated line profiles. 

The main conclusions of our study are the following: 

- For stars with significant temperature inhomogeneities a re- 
liable magnetic field reconstruction with ZDI is essentially 
impossible when temperature is not taken into account dur- 
ing mapping of magnetic spots. If a multi-component nature 
of the active-star surface is not taken into account in mag- 
netic mapping, the magnetic field strength inside cool spots 
is grossly underestimated and the overall geometry of the 
field is not recovered correctly. 

- Our tests also showed that even if the temperature distribu- 
tion is homogeneous there will still be a discrepancy between 
reconstructed magnetic field maps depending on if Stokes 
/ is used in mapping or not. The temperature reconstruc- 
tion itself was accurate to within 100 K and should therefore 
not have caused the relatively large difference in the recon- 
structed magnetic field maps. Instead, we found that addi- 
tional information contained in the Zeeman broadening and 
intensification of Stokes / profiles helps to constrain mag- 
netic maps. 

- The shape of Stokes / line profiles is affected by a magnetic 
field making them broader and rotationally modulated even 
in the absence of temperature inhomogeneities. It may there- 
fore be important to take a magnetic field into account dur- 
ing temperature mapping of active cool stars with especially 
strong fields. 

- According to our results, simultaneous recovery of magnetic 
field and temperature is feasible with a self-consistent inter- 
pretation of the Stokes / and V profiles of individual spectral 
lines. We did not find significant crosstalks between temper- 
ature and magnetic maps. Instead, reconstruction of mag- 
netic field is greatly improved when appropriate local tem- 
perature variations are taken into account. Thus, there are 
no reasons whatsoever to keep temperature and magnetic in- 



versions separate, as done by the majority of previous and 
current ZDI studies. 

- Magnetic mapping using Stokes / and V suffers from a num- 
ber of systematic artifacts. For example, it is difficult to dis- 
tinguish radial and meridional fields at low latitudes and the 
strength of the reconstructed radial field generally decreases 
with decreasing latitude. At the same time, azimuthal field 
is reliably recovered for a wide range of spot latitudes and 
does not suffer from a crosstalk with either radial or merid- 
ional field. 

- When the linear polarization data are incorporated in the ZDI 
inversions, the reconstructed field topology becomes notice- 
ably more accurate compared to the inversions with partial 
Stokes data sets. In particular, by including all four Stokes 
parameters one can avoid crosstalks between different field 
components. 
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